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Abstract 

Early photoionization of the intergalactic medium is discussed in 
a nearly model-independent way, in order to investigate whether early 
structures corresponding to rare Gaussian peaks in a CDM model can 
photoionize the intergalactic medium sufficiently early to apprecia- 
bly smooth out the microwave background fluctuations. We conclude 
that this is indeed possible for a broad range of CDM normalizations 
and is almost inevitable for unbiased CDM, provided that the bulk of 
these early structures are quite small, no more massive than about 
IQ^Mq. Typical parameter values predict that reionization occurs 
around z = 50, thereby suppressing fluctuations on degree scales while 
leaving the larger angular scales probed by COBE reasonably unaf- 
fected. However, for non-standard CDM, incorporating mixed dark 
matter, vacuum density or a tilted primordial power spectrum, early 
reionization plays no significant role. 
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1 Introduction 



The first quantitative predictions of cosmic microwave background anisotropics 
in cold dark matter (CDM)-dominated cosmological models recognized that 
reionization by rare, early-forming objects could play a role in suppressing 
temperature fluctuations on small angular scales (Bond & Efstathiou 1984; 
Vittorio k Silk 1984). Now that the COBE DMR experiment has detected 
fluctuations on large angular scales (Smoot et al. 1992) at a level (within a 
factor of two) comparable to that predicted by CDM models, it is especially 
relevant to examine whether reionization can affect the degree scale searches 
that are currently underway. 

Cold dark matter models are generally characterized by a late epoch of 
galaxy formation. However, the smallest and oldest objects first go nonlinear 
at relatively large redshift. In this paper we investigate, for a wide range 
of CDM normalizations, power spectra and efficiency parameters, whether 
reionization associated with energy injection by early forming dwarf galaxies 
can reionize the universe sufficiently early to smooth out primordial CBR 
temperature fluctuations. 

Although we go into some detail in the appendix to make estimates 
of a certain efficiency parameter, our overall treatment is fairly model- 
independent, and can be used as a framework within which to compare 
various photoionization scenarios. Our basic picture is roughly the follow- 
ing: An ever larger fraction of the baryons in the universe falls into non- 
linear structures and forms galaxies. A certain fraction of these baryons 
form stars or quasars which emit ultraviolet radiation, and some of this ra- 
diation escapes into the ambient intergalactic medium (IGM) and ends up 
photoionizing and heating it. Due to cooling losses and recombinations, the 
net number of ionizations per UV photon is generally less than unity. 

Apart from photoionization, early galaxies can also ionize the IGM through 
supernova driven winds, an ionization mechanism that will not be treated in 
this paper. Although such winds can ionize the IGM by z = 5, early enough 
to explain the absence of a Gunn- Peterson effect (Tegmark et al. 1993), the 
relatively low velocities of such winds makes them unable to distribute the 
released energy throughout space at redshifts early enough (by z ^ 50) to 
measurably affect the CBR. 

Our approach will be to first write the ionization fraction of the IGM as 
a product of a number of factors, and then discuss the value of each of these 
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factors in more detail. Let us write 



X — fsfuvppfi 



uvpp J loni 



where 

X 

fs 

fuvpp 
fion 



fraction of IGM that is ionized, 

fraction of baryons in nonUnear structures, 

UV photons emitted into IGM per proton in nonhnear structures, 
net ionizations per emitted UV photon. 



Let us first consider the case where the UV photons are produced by stars, 
and return to the quasar case later. Using the fact that a fraction 0.0073 of 
the rest mass is released in stellar burning of hydrogen to helium, we obtain 



uvpp 



0.0073 



13.6eV 



Ih fburn fuv fe 



where 



Jh = mass fraction hydrogen in IGM, 

fburn = niass fraction of hydrogen burnt, 

fuv = fraction of energy released as UV photons, 

fesc = fraction of UV photons that escape from galaxy. 

We will take the primordial mass fraction of helium to be 24%, i.e. fn 
76%. Now define the net efficiency 



fnet — fburn fuv fesc fu 



and equation (Q) becomes 



X ~ 3.8 X 10^ fs 



(3) 



The key feature to note about this expression is that since 3.8 x 10^ is such 
a large number, quite modest efficiencies fnet still allow x to become of or- 
der unity as soon as a very small fraction of the baryons are in galaxies. 
As will be seen in the next section, this means that reionization is possible 
even at redshifts far out in the Gaussian tail of the distribution of forma- 
tion redshifts, at epochs long before those when the bulk of the baryons go 
nonlinear. This appears to have been first pointed out by Couchman and 
Rees (1986). 
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2 The Mass Fraction in Galaxies 



In this section, we will discuss the parameter fg. Assuming the standard 
PS theory of structure formation (Press & Schechter 1974), the fraction of 
all mass that has formed gravitationally bound objects of total (baryonic 
and non-baryonic) mass greater than M at redshift z is the integral of the 
Gaussian tail, 

Sc 



fs = erfc 



V2a{M, 



(4) 



where the complementary error function erfc(a;) = 2tt^^^'^ e'^^ du and 
a{M, z) is the r.m.s. mass fluctuation in a sphere containing an expected 
mass M at redshift z. is given by top-hat filtering of the power spectrum 
as 

" sin kvQ cos kr^ ^ 



a{M,zf oc / P{k) 
Jo 



(kro)^ {kroY 



dk, (5) 



where P{k) is the power spectrum at redshift z and ro is given by ^t^VqP = 
M, p = being the density of the universe at redshift z. Although 

this approach has been criticized as too simplistic, numerical simulations 
(Efstathiou et al. 1988; Efstathiou & Rees 1988; Carlberg & Couchman 
1989) have shown that it describes the mass distribution of newly formed 
structures remarkably well. Making the standard assumption of a Gaussian 
density field, Blanchard et al. (1992) have argued that it is an accurate 
description at least in the low mass limit. Since wc arc mainly interested in 
extremely low masses such as IO^Mq, it appears to suffice for our purposes. 

For our middle-of-the road estimate, we choose 5c = 1.69, which is the 
linearly extrapolated overdensity at which a spherically symmetric pertur- 
bation has collapsed into a virializcd object (Gott & Rees 1975; Efstathiou 
et al. 1988; Brainerd & Villumsen 1992). We take 6c = 1.44 (Carlberg & 
Couchman 1989) for the optimistic estimate, although the even lower value 
5c = 1.33 has been discussed (Efstathiou Sz Rees 1988), and 6c = 2.00 (Gelb 
&: Bcrtschinger 1992) for the pessimistic estimate. (Here and throughout 
this paper, parameter choices are referred to as optimistic if they permit 
earlier reionization.) 

The fact that a{M,z) — oo as M — implies that if we consider 
arbitrarily small scales, then all dark matter is in non-linear structures. 
Thus if no forces other than gravity were at work, so that the baryons 
always followed the dark matter, we would simply have = 1 at all z. 
However, it is commonly believed that galaxies correspond only to objects 
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that are able to cool (and fragment into stars) in a dynamical time or a 
Hubble time (Binney 1977; Rees & Ostriker 1977; Silk 1977; White & Rees 
1978). The former applies to ellipticals and bulges, the latter to disks. Let 
us define the virialization redshift (1 + Zvir) = {V^/Sc)(^{Mc,0), where Mc 
is some characteristic cutoff mass which is the total mass (baryonic and 
dark) of the first galaxies to form. Zmr is roughly the redshift at which the 
bulk of all baryons goes non-linear. Using equation and the fact that 
cr(M, z) = a{M, 0)/(l + z) in the linear regime of CDM, we thus have 



fs = erfc 



1 + z 



1+ Zy 



(6) 



A common assumption is that ~ 10^ Mq, roughly the Jeans mass at re- 
combination. Blanchard et al. (1992) examine the interplay between cooling 
and gravitational collapse in considerable detail, and conclude that the first 
galaxies to form have masses in the range 1O^M0 to 10® Mq, their redshift 
distribution still being given by equation (|6|), whereas Couchman &; Rees 
(1986) argue that the first galaxies to form may have had masses as low as 
IO^Mq. 

As our CDM power spectrum today, we will use that given by BBKS 
(Bardeen et al. 1986) and an n = 1 Harrison-Zel'dovich primordial spectrum: 

^ ( g"Mn(l + 2.34g) V 

V[l + 3.89g + (16.1g)2 + (5.46g)3 + (6.71g)4]i/V ' 

where q = k /[h?VlQ\sA.^c~^]. Throughout this paper, we will take JIq = 1- 
Evaluating the u^-integral in equation (^) numerically yields 

t7(lO^M0,O) « 33.76"^ 

for h = 0.8 and 

£7(10® M0,O) f« 13.66"^ 

for h = 0.5, where the so called bias factor h = a{8h~^Mpc,0)~^ has been 
estimated to lie between 0.8 (Smoot et al. 1992) and 2.5 (Bardeen et al. 
1986). Our pessimistic, middle-of-the-road and optimistic CDM estimates 
of Zmr are given in Table ||, and the dependence of z^r on is plotted 
in Figure |l[ This figure also shows three alternative models of structure 
formation: CDM with cosmological constant (Efstathiou et al. 1992); tilted 
CDM (Cen et al. 1992) and MDM, mixed hot and dark matter (Shafi & 
Stecker 1984; Schaefer & Shafi 1992; Davis et al. 1992; Klypin et al. 1993). 
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For the model with cosmological constant, we have taken a flat universe with 
h = 0.5, = 0.4 and Aq = 0.6. For the tilted model, the power spectrum 
P{k) is simply multiplied by a factor A;"~^, where we have taken n = 0.7. 
For the tilted case, equation (|6|) still applies. For the MDM case, however, 
perturbations in the cold component grow slower than linearly with the scale 
factor (1 + z)~^ and equation @ is not valid. For the low masses we are 
considering, we have (Bond &: Szalay 1983) 



(tmdm{Mc,z) 



(Tmdm{Mc,0) 



where 



""4 



v/25 - 2AnHDM - 1 



Using the parameters from Davis et al. (1992), who take $1 = 1 and ^hdm 
0.3, the MDM version of equation dq) becomes 



fs = erfc 



1 + z 

1 ~\~ ^vii 



where a ~ 0.8 and we redefine z^r by 



(l + ^^, 



^a(M„0) 



-1 1/q 



with (t(Mc, 0) referring to a pure CDM power spectrum. 

For the A case, perturbations grow approximately linearly until the uni- 
verse becomes vacuum dominated at z ~ fig ^ — 1 = 1.5, after which their 
growth slowly grinds to a halt. A numerical integration of the Friedmann 
equation and the equation for perturbation growth using h = 0.5, VLq = 0.4 
and Ao = 0.6 gives 

,^(Me,0) 



oa{Mc,z) 



1.2- 



(1 + ^) 



for z 3> 3. 

Since our A-model yields a value of Zm^ 
we will omit the former from future plots. 



very similar to our tilted model. 



3 Efficiency Parameters 

In this section, we will discuss the various parameters that give fnet when 
multiplied together. The conclusions are summarized in Table |2|. 
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Mixed 


Tilted 


Lambda 


Pess. 


Mid. 


Opt. 


Me 


W^Mq 




IO'^Mq 


1O«M0 


W^Mq 


IO^Mq 


Model 


MDM 


Tilted 


Lambda 


CDM 


CDM 


CDM 


h 


0.5 


0.5 


0.5 


0.5 


0.5 


0.8 


b 


1 


1 


1 


2 


1 


0.8 




1.69 


1.69 


1.69 


2.00 


1.69 


1.44 




2.9 


10.1 


8.4 


4.8 


17.2 


41.4 



Table 1: Galaxy formation assumptions 





Pess. 


Mid. 


Opt. 


fburn 


0.2% 


1% 


25% 


fesc 


10% 


20% 


50% 


fuv 


5% 


25% 


50% 


fion 


10% 


40% 


95% 


fnet 


1 X lO-'* 


2 X 10-* 


6 X 10-^ 


fuvpp 


4 


190 


24,000 



Table 2: Efficiency parameters used 
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UV source 


Spectrum P(z^) 


fuv 


(Euv) 




0"18 


03 star 


T = 


50, OOOK Planck 


0.57 


17.3 eV 


28,300K 


2.9 


06 star 


T = 


40, OOOK Planck 


0.41 


16.6 eV 


23,400K 


3.4 


09 star 


T = 


30, OOOK Planck 


0.21 


15.9 eV 


18,000K 


3.9 


Pop. Ill star 


T = 


50, OOOK Vacca 


0.56 


18.4 eV 


36,900K 


2.2 


Black hole, QSO 


a = 


1 power law 




18.4 eV 


37,400K 


1.7 


? 


a = 


2 power law 




17.2 eV 


27,800K 


2.7 


? 


a = 


power law 




20.9 eV 


56,300K 


0.6 


? 


T = 


100, OOOK Planck 


0.89 


19.9 eV 


49,000K 


1.6 



Table 3: Spectral parameters 



fburm the fraction of galactic hydrogen that is burnt into helium dur- 
ing the early life of the galaxy (within a small fraction of a Hubble after 
formation), is essentially the galactic metallicity after the first wave of star 
formation. Thus it is the product of the fraction of the hydrogen that forms 
stars and the average metallicity per star (weighted by mass). This depends 
on the stellar mass function, the galactic star formation rate and the final 
metallicities of the high-mass stars. For our middle-of-the-road estimate, we 
follow Miralda-Escude &: Ostriker (1990) in taking ftum = 1%, half the solar 
value. An upper limit to fbum is obtained from the extreme scenario where 
all the baryons in the galaxy form very massive and short-lived stars with 
M « 30Mq , whose metallicity could get as high as 25% (Woosley & Weaver 
1986). Although perhaps unrealistic, this is not ruled out by the apparent 
absence of stars with such metallicities today, since stars that massive would 
be expected to collapse into black holes. 

In estimating f^sc, the fraction of the UV photons that despite gas and 
dust manage to escape from the galaxy where they are created, we follow 
Miralda-Escude & Ostriker (1990). 

For fuy, the fraction of the released energy that is radiated above the 
Lyman limit, we also follow Miralda-Escude & Ostriker (1990). The upper 
limit refers to the extreme 3OM0 scenario mentioned above. For reference, 
the values of fuv for stars with various spectra are given in Table |3|, together 
with some other spectral parameters that will be defined and used in the 
appendix. All these parameters involve spectral integrals, and have been 
computed numerically. 

The parameter fion is estimated in the appendix. 
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An altogether different mechanism for converting the baryons in non- 
linear structures into ultraviolet photons is black hole accretion. If this 
mechanism is the dominant one, equation (|2|) should be replaced by 

fuvpp ~ ^ 13 6eV^ fbhfaccfuvfesc: 

where 

fljfi = mass fraction of nonlinear structures that end up as black holes, 
face = fraction of rest energy radiated away during accretion process, 
fuv = fraction of energy released as UV photons, 
fese = fraction of UV photons that escape from host galaxy. 

There is obviously a huge uncertainty in the factor However, the ab- 
sence of the factor 0.0073 x 0.76 compared to equation (||) means that the 
conversion of matter into radiation is so much more efficient that the black 
hole contribution might be important even if f^h is quite small. For instance, 
faec = 10% and fese = W0% gives fuvpp ~ lO^fbhfuv, which could easily ex- 
ceed the optimistic value fuvpp ~ 24, 000 for the stellar burning mechanism 
in Table |. 

In Figure the ionization fraction xi^) is plotted for various parameter 
values using equations (|l|) and @. It is seen that the ionization grows quite 
abruptly, so that we may speak of a fairly well-defined ionization redshift. 
Let us define the ionization redshift Zion as the redshift when x becomes 0.5, 
i.e. 

1 + Zion = (1 + ^;i,ir)erfc"^ — — . (7) 

y ^ J uvpp J ion I 

This dependence of Zion of the efficiency is shown in Figure ^ for our var- 
ious galaxy formation scenarios. It is seen that the ionization redshift is 
fairly insensitive to the net efficiency, with the dependence being roughly 
logarithmic for fnet > 0.0001. 



4 Scattering History 

For a given ionization history xi^), the Thomson opacity out to a redshift 
z, the probability that a CBR photon is Thomson scattered at least once 
after z, is 

Ps{z) = l-e-^^'\ 
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where the optical depth for Thomson scattering is given by 



T[Z 



Stt 



1 



Hpcat 



0.057/iQb, 



where we have taken the mass fraction of hehum to be fne ~ 24% and 
assumed ~ x, i.e. that hehum never becomes doubly ionized and that 
the fraction that is singly ionized equals the fraction of hydrogen that is 
ionized. The latter is a very crude approximation, but makes a difference 
of only 6%. We assume that = 1 throughout this paper, ilf, denotes 
the density of the intergalactic medium divided by the critical density, and 
is usually assumed to equal Of,, the corresponding density of baryons. the 
probability that a CBR photon is Thomson scattered at least once after the 
standard recombination epoch at z ~ 10^. 

The profile of the last scattering surface is given by the so called visibility 
function 

az 

which is the probability distribution for the redshift at which a photon last 
scattered. An illuminating special case is that of complete ionization at all 
times, i.e. x(-z) = 1, which yields 



Ps{z) = 1 - exp 



--T* 

3 



(1 + ^) 



3/2 



1 — exp 



92 



3/2' 



(8) 



for z ^ I and hQ^ = 0.03. Hence we see that in order for any significant 
fraction of the CBR to have been rescattered by reionization, the reionization 
must have occurred quite early. Figures ^ and ^ show the opacity and last- 
scattering profile for three different choices of hQf,. In the optimistic case 
/if2j, = 0.1, it is seen that even as low an ionization redshift as z^Qn = 30 
would give a total opacity Pg ~ 50%. In Figures ^ and |^, we have replaced 
z by the angle subtended by the horizon radius at that redshift. 



9{z) 



2 arctan 



2(vrTi- 1) 



which is the largest angular scale on which Thomson scattering at z would 
affect the microwave background radiation. In Figure 0, we have plotted the 
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angular visibility function dPs/ d{—9) instead of dPs/dz, so that the curves 
are probability distributions over angle instead of redshift. 

In the sudden approximation, the ionization history is a step function 

Xiz) = 9{zion - z) 

for some constant zion-, and as was discussed in Section this models the 
actual ionization history fairly well. In this approximation, the visibility 
functions are identical to those in Figures ^ and ^ for z < Zion, but vanish 
between Zion and the recombination epoch at 2: ~ 10^. Figure H], which is 
in a sense the most important plot in this paper, shows the total opacity 
Psi,Zion^ as a function of fnet for a variety of parameter values, as obtained 
by substituting equation (0) into (|^). As can be seen, the resulting opacity 
is relatively insensitive to the poorly known parameter fnet, and depends 
mainly on the structure formation model (i.e. z^r) and the cosmological 
parameter hil.f). 

Thomson scattering between CBR photons and free electrons affects not 
only the spatial but also the spectral properties of the CBR. It has long 
been known that hot ionized IGM causes spectral distortions to the CBR, 
known as the Sunyaev-Zel'dovich effect. A useful measure of this distortion is 
the Comptonization y-parameter (Kompaneets 1957; Zel'dovich & Sunyaev 
1969; Stebbins & Silk 1986; Bartlett & Stebbins 1991) 

yc= [ f — — ^) neatc dt, 

where the integral is to be taken from the reionization epoch to today. Let 
us estimate this integral by making the approximation that the IGM is cold 
and neutral until a redshift Zion, at which it suddenly becomes ionized, and 
after which it is remains completely ionized with a constant temperature T. 
Then for = 1, Zion ^ 1, we obtain 

where = T/IO^K and rieo, the electron density today, has been computed 
as before assuming a helium mass fraction of 24% that is singly ionized. 
Substituting the most recent observational constraint from the COBE FI- 
RAS experiment, yc < 2.5 x 10^^ (Mather et al. 1994), into this expression 
yields 

Znon < 554i. , 

so all our scenarios are consistent with this spectral constraint. 
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5 Discussion 



A detailed discussion of how reionization affects the microwave background 
anisotropics would be beyond the scope of this paper, so we will merely re- 
view the main features. If the microwave background photons are rescattered 
at a redshift z, then the fluctuations we observe today will be suppressed 
on angular scales smaller than the angle subtended by the horizon at that 
redshift. This effect is seen in numerical integrations of the linearized Boltz- 
mann equation {e.g. Bond &: Efstathiou 1984; Vittorio & Silk 1984), and 
can be simply understood in purely geometrical terms. Suppose we detect 
a microwave photon arriving from some direction in space. Where was it 
just after recombination? In the absence of reionization, it would have been 
precisely where it appears to be coming from, say 3000 Mpc away. If the 
IGM was reionized, however, the photon might have originated somewhere 
else, scattered off of a free electron and then started propagating towards us, 
so at recombination it might even have been right here. Thus to obtain the 
observed anisotropy, we have to convolve the anisotropics at last scattering 
with a window function that incorporates this smoothing effect. Typical 
widths for the window function appropriate to the last scattering surface 
range from a few arc-minutes with standard recombination to the value of 
a few degrees that we have derived here for early reionization models. 

In addition to this suppression on sub-horizon scales, new fluctuations 
will be generated by the first order Doppler effect and by the Vishniac 
effect. The latter dominates on small angular scales and is not included 
in the linearized Boltzmann treatment because it is a second order effect. 
The current upper limit on CBR fluctuations on the 1 arcminute scale of 
AT/T < 9 X 10~^ (Siibrahmanyan et al. 1993) provides an interesting con- 
straint on reionization histories through the Vishniac effect. In fact, accord- 
ing to the original calculations (Vishniac 1987), this would rule out most of 
the reionization histories in this paper. However, a more careful treatment 
(Hu et al. 1994) predicts a Vishniac effect a factor of five smaller on this 
angular scale, so all reionization histories in this paper are still permitted. 

The COBE DMR detection of AT/T has provided a normalization for 
predicting CBR anisotropies on degree scales. Several experiments are un- 
derway to measure such anisotropies, and early results that report possible 
detections have recently become available from experiments at the South 
Pole (Meinhold &: Lubin 1991; Shuster et al. 1993) and at balloon altitudes 
(Devlin et al. 1992; Meinhold et al. 1993; Shuster et al. 1993). There is some 
reason to believe that these detected signals are contaminated by galactic 
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emission. Were this the case, the inferred CBR upper hmits to fluctuations 
on degree scales might be lower than those predicted from COBE extrap- 
olations that adopt the scale-invariant power spectrum that is consistent 
with the DMR result and is generally believed to be the most appropriate 
choice on large scales from theoretical considerations {e.g. Gorski et al. 1993; 
Kashlinsky 1992). In the absence of such contaminations, the detected fluc- 
tuations in at least some degree-scale experiments are, however, consistent 
with the COBE extrapolation {e.g. Jubas & Dodelson 1993). The varia- 
tion from field to field, repeated on degree scales, also may argue either 
for galactic contamination or else for unknown experimental systematics, 
or even non-Gaussian fluctuations. The results of other recent experiments 
such as ARGO (de Bernardis et al. 1993), PYTHON (Dragovan et al. 1993) 
and MSAM (Cheng et al. 1993) have reinforced the impression that the 
experimental data is not entirely self-consistent, and that some form of sys- 
tematic errors may be important. 

The controversy over the interpretation of the degree-scale CBR fluctua- 
tions makes our reanalysis of the last scattering surface particularly timely. 
We have found that canonical dark matter, tailored to provide the 10 de- 
gree CBR fluctuations detected by the COBE DMR experiment, results in 
sufficiently early reionization (before z ~ 50) over a fairly wide range of 
parameter space, to smooth out primordial degree-scale flTictuations. Our 
middle-of-the-road model produces suppression by roughly a factor of two; 
it is difficult, although not impossible, to obtain a much larger suppres- 
sion. This smoothing, because it is of order unity in scattering optical 
depth, is necessarily inhomogeneous. We predict the presence of regions 
with large fluctuations and many "hot spots" and "cold spots" , correspond- 
ing to "holes" in the last-scattering surface, as well as regions with little 
small-scale power where the last scattering is more efficient. The detailed 
structure of the CBR sky in models with reionization will be left for future 
studies. Here we simply conclude by emphasizing that anomalously low val- 
ues of AT/T over degree scales are a natural corollary of reionization at 
high redshift. 



A Appendix: The Efficiency Parameter fion 

In this appendix, we will discuss the parameter fion, and see that it rarely 
drops below 30%. We will first discuss the thermal evolution of intergalactic 
hydrogen exposed to a strong UV fiux, and then use the results to write down 
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a differential equation for the volume fraction of the universe that is ionized, 
subject to point sources of UV radiation that switch on at different times. 
Wc will see that photoionization is so efficient within the ionized regions of 
the IGM that quite a simple equation can be given for the expansion of the 
ionized regions. 

The evolution of IGM exposed to ionizing radiation has been discussed 
by many authors. Important early work includes that of Arons & McCray 
(1970), Bergeron & Salpeter (1970) and Arons & Wingert (1972). The main 
novelty of the treatment that follows is that whereas previous treatments 
focus on late {z < 5) epochs, when various simplifying approximations can 
be made because the recombination and Compton rates are low, we are 
mainly interested in the case 50 < z < 150. We show that IGM exposed to 
a strong UV flux rapidly approaches a quasistatic equilibrium state, where 
it is almost fully ionized and the temperature is such that photoionization 
heating exactly balances Compton cooling. This simplifies the calculations 
dramatically, since the entire thermal history of the IGM can be summarized 
by a single function xi^)^ the volume fraction that is ionized. Thus a fraction 
x{z) is ionized and hot (with a temperature that depends only on z, not on 
when it became ionized), and a fraction 1 — x{z) is neutral and cold. 

In the first section, we justify this approximation. In the second section, 
we derive a differential equation for the time-evolution of x as well as a 
useful analytic estimate of fion- 

A.l Intergalactic Stromgren Spheres 

Let X denote the ionization fraction in a small, homogeneous volume of in- 
tergalactic hydrogen, i.e. x = nnii / {nni + nnii)- {x is not to be confused 
with X, the volume fraction in ionized bubbles.) When this IGM is at tem- 
perature r, exposed to a density of rj UV photons per proton, the ionization 
fraction x evolves as follows: 



where Hq^{1 + z) ^ times the rates per baryon for photoionization, coUi- 
sional ionization and recombination are given by 



dx 1 + z 



\pi(l - x) + Xcix{l - x) - A^g],x^ 



(9) 




Xt-ec ~ Q.717hnbT^ ' 1.808 - 0.51nr4 + 0.187r4 



(10) 
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and T4 = r/lO'^K. ais is the thermally averaged photoionization cross 
section in units of 10~^^cm^, and has been computed in Table ^ for various 
spectra using the differential cross section from Osterbrock (1974). The 
collisional ionization rate is from Black (1981). The recombination rate is 
the total rate to all hydrogenic levels (Seaton 1959). 

Below we will see that in the ionized Stromgren bubbles that will appear 
around the galaxies or quasars, the photoionization rate is so much greater 
than the other rates that to a good approximation, equation (^) can be 
replaced by the following simple model for the IGM: 

• It is completely ionized (x = 1). 

• When a neutral hydrogen atom is formed through recombination, it is 
instantly photoionized again. 

Thus the only unknown parameter is the IGM temperature, which de- 
termines the recombination rate, which in turn equals the photoionization 
rate and thus determines the rate of heating. 

Let us investigate when this model is valid. Near the perimeter of an 
ionized Stromgren sphere of radius r surrounding a galaxy, the number of 
UV photons per proton is roughly 

V 



where Suv is the rate at which UV photons leave the galaxy. For an 05 
star, the photon flux above the Lyman limit is approximately 3.1 x lO^^s"^ 
(Spitzer 1968), so if each solar masses of baryons in a galaxy leads to 
production of a UV flux equivalent to that of an 05 star, then 

n, r^tj fescMg n 1\ 

/i^rf iV(l + z)3 ' ^ ^ 

inside the sphere, where ri = r/lMpc and Mq = M/10^Mq. When a 
fraction fg of all matter has formed galaxies of a typical total (baryonic 
and dark) mass M, then in the absence of strong clustering, the typical 
separation between two galaxies is 



R 



Af /i5kpc\ / Mq 



fspj \l + Z 



Km 



where Mq = M/W^Mq. Thus r continues to increase until r ^ R, and 
spheres from neighboring galaxies begin to overlap. We are interested in 
the regime where z < 150. Substituting this and equation (^) into (p^, 
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we see that Xpi » Xd and Xpi » Xrec for any reasonable parameter values. 
Hence we can neglect collisional ionization in equation (|9|). Since Xpi ^ 1, 
the photoionization timescale is much shorter than the Hubble timescale, so 
equation (^) will quickly approach a quasistatic equilibrium solution where 
the recombination rate equals the photoionization rate, i.e. 

-, Xrec -, 
X PS 1 1. 

In conclusion, the simple x = 1 model is valid for all parameter values in 
our regime of interest. 

When a hydrogen atom gets ionized, the photoelectron acquires an aver- 
age kinetic energy of |/cT*, where T* is defined by ^kT* = {E^v) — 13.6eV, 
and {Euv) is the average energy of the ionizing photons (see Table ^). 

Since the timescale for Coulomb collisions is much shorter than any other 
timescales involved, the electrons and protons rapidly thermalize, and we can 
always assume that their velocity distribution is Maxwellian, corresponding 
to some well-defined temperature T. Thus shortly after the hydrogen gets 
photoionized, after the electrons have transferred half of their energy to the 
protons, the plasma temperature is T = \T* . 

The net effect of a recombination and subsequent photoionization is to 
remove the kinetic energy of the captured electron, say ^kTr]rec{T), from the 
gas and replace it with |A;T*, the kinetic energy of the new photoelectron. 
Since the recombination cross section is approximately proportional to v~'^, 
slower electrons are more likely to get captured. Hence the mean energy of 
the captured electrons is slightly lower than ^kT, i.e. r]rec{T) is slightly less 
than unity (Osterbrock 1974). We compute rj^edT) using Seaton (1959). 
The complication that rfreciT) ^ 1 turns out to be of only marginal im- 
portance: r/rec(10'*K) ~ 0.8, which only raises the equilibrium temperatures 
calculated below by a few percent. 

The higher the recombination rate, the faster this effect will tend to 
push the temperature up towards T* . The two dominant cooling effects 
are Compton drag against the microwave background photons and cooling 
due to the adiabatic expansion of the universe. Line cooling from collisional 
excitations can be neglected, since the neutral fraction 1 — x ~ 0. Combin- 
ing these effects, we obtain the evolution equation for the IGM inside of a 
Stromgren bubble: 



dT 2 l + z 



d{-z) l + z ^/TTThz 



^^comp {Tebr - T) + -XreciT)[T,br " r]rec{T)T] 

(12) 
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where 

= 1^ hi?) id-J' + ^) =^ 0.00418/.-(l + .) 

is (1 + z)^^ times the Compton cooUng rate per Hubble time and Tcbr = 
TcbrfiC^ + z). We have taken T^brfi ~ 2.726K (Mather et al. 1994). The factor 
of ^ in front of the Xrec term is due to the fact that the photoelectrons share 
their acquired energy with the protons. The average energy of the ionizing 
photons is given by the spectrum P{i^) as (Euv) = h{v), where 

Here the photoionization cross section (t(i^) is given by Osterbrock (1974). 
Note that, in contrast to certain nebula calculations where all photons get 
absorbed sooner or later, the spectrum should be weighted by the photoion- 
ization cross section. This is because most photons never get absorbed in 
the Stromgren regions (only in the transition layer), and all that is rele- 
vant is the energy distribution of those photons that do. P{v) is the energy 
distribution (W/Hz), not the number distribution which is proportional to 
P{v)/v. 

The spectral parameters {E^v) and T* are given in Table |3| for some 
selected spectra. A Planck spectrum P{v) oc / (^ef^^/^'^ — \j gives quite a 
good prediction of T* for stars with surface temperatures below 30, OOOK. 
For very hot stars, more realistic spectra (Vacca, 1993) have a sharp break 
at the Lyman limit, and fall off much slower above it, thus giving higher 
values of T* . As seen in Table ^ an extremely metal poor star of surface 
temperature 50, OOOK gives roughly the same T* as QSO radiation. The 
only stars that are likely to be relevant to early photoionization scenarios 
are hot and short-lived ones, since the universe is only about 10'' years old 
at z = 100, and fainter stars would be unable to inject enough energy in so 
short a time. Conceivably, less massive stars could play a dominant role later 
on, thus lowering T* . However, since they radiate such a small fraction of 
their energy above the Lyman limit, very large numbers would be needed, 
which could be difficult to reconcile with the absence of observations of 
Population HI stars today. If black holes are the dominant UV source, the 
stellar spectra of Table ^ are obviously irrelevant. A power law spectrum 
P(z/) oc with a = 1 fits observed QSO spectra rather well in the vicinity 
of the Lyman limit (Cheney & Rowan-Robinson 1981; O'Brien et al. 1988), 
and is also consistent with the standard model for black hole accretion. 
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Numerical solutions to equation (|T2|) are shown in Figure and it is 
seen that the temperature evolution separates into three distinct phases. 
In the first phase, the IGM is outside of the Stromgren regions, unexposed 
to UV radiation, and remains cold and neutral. In the second phase, the 
IGM suddenly becomes ionized, and its temperature instantly rises to ^T*. 
After this, Compton cooling rapidly reduces the temperature to a quasi- 
equilibrium value of a few thousand K. After this, in the third phase, T 
changes only quite slowly, and is approximately given by setting the expres- 
sion in square brackets in equation (12) equal to zero. This quasi-equilibrium 
temperature is typically many times lower than T* , since Compton cooling 
is so efficient at the high redshifts involved. 



A. 2 The Expansion of Stromgren Regions 

This rapid approach to quasi-equilibrium, where the IGM "loses its memory" 
of how long ago it became part of a Stromgren region, enables us to construct 
a very simple model for the ionization history of the universe. At redshift z, 
a volume fraction xi^) of the universe is completely ionized and typically has 
a temperature of a few thousand K. The ionized part need not consist of non- 
overlapping spheres; it can have any topology whatsoever. The remainder 
is cold and neutral. 

Between the ionized and neutral regions is a relatively thin transition 
layer, where the IGM becomes photoionized and its temperature adjusts to 
the quasistatic value as in Figure As this IGM becomes part of the hot 
and ionized volume, the transition layer moves, and x{^) increases^. 

As long as X < 1; ^-ll UV photons produced are absorbed instantly to a 
good approximation. Thus the rate at which UV photons are released is the 
sum of the rate at which they are used to counterbalance recombinations 
inside the hot bubbles and the rate at which they are used to break new 
ground, to increase x- Thus 

^""PP d t ^ ^ dt ' 

^We are tacitly assuming that the UV luminosity of the galaxy that creates each 
Stromgren sphere never decreases. Although obviously untrue, this is in fact an excellent 
approximation, since these early dwarf-galaxies correspond to perturbations far out in the 
Gaussian tail. Since fs{z) grows so dramatically as the redshift decreases and we move 
from five sigma to four sigma to three sigma, etc., almost all galaxies in existence at a 
given redshift are in fact very young, so that older ones that have begun to dim can be 
safely neglected. 
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where a^^^ (T) is the total recombination rate to all hydrogenic levels except 
the ground state^. Changing the independent variable to redshift and using 
equation (0), we find that this becomes 



+ A(2) -i±^;^ = A f J!fM^ ) exp 



d{-z) ^'^Vl + f^o^ V^\l + z^, 
Here 



1 + z 

1 H~ ^vir 



(13) 



1.04-0.51nr4 + 0.19r2/^ 



is Hq ^(1 + z) ^ times the total recombination rate per baryon to all hy- 
drogenic levels except the ground state. The fit is to the data of Spitzer 
(1968) and is accurate to within 2% for 30K < T < 64, OOOK. A^ec is to be 
evaluated at the quasi-equilibrium temperature T{z) discussed above. 

Using the values in Table ^ for the pessimistic, middle-of-the-road and 
optimistic estimates, the parameter fuvpp equals roughly 4, 190 and 24,000, 
respectively. 

In the absence of photon waste through recombination, equation equa- 
tion (p^) would have the solution x*iz) = fuvppfsiz), so the ionization effi- 
ciency is 

fioniz) = X{z)/X*{z). 

Since equation (|l^) is linear in x the initial data is x = ^-t some 
redshift, it is readily seen that the solution xiz) is proportional to fuvpp, the 
constant in front of the source term. Combining these last two observations, 
we see that fion is independent of fuvpp and hence independent of the poorly 
known parameters fbum, fuv and fesc- 

Plots of fion{z) from numerical solutions of equation (|T3| ) are shown in 
Figure |^ for various parameter values, and it is seen that the dependence 
on z is generally quite weak. Let us make use of this fact by substituting the 
Ansatzxiz) = f ion{z) fuvpp fs{z) into equation (||), and setting /•^^(z) 0. 
Using equation ^ and an asymptotic approximation for the error function, 

''a reionization directly to the ground state produces a UV photon that usually propa- 
gates uninterrupted through the highly ionized Stromgren region, and then ionizes another 
atom in the transition layer between the expanding Stromgren region and its cold and neu- 
tral surrounding. Thus recombinations directly to the ground state were included in the 
above calculation of the quasi-equilibrium temperature of the Stromgren bubbles, since 
the resulting UV photons could be considered lost from the latter. Here, on the contrary, 
recombinations directly to the ground state should not be included, since the UV photons 
they produce are not wasted from an energetics point of view. 
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we obtain ^ 

1 + o.48AreUi + z^ir)yvT+n^ 

independent of fuvpp-i which agrees to within 10% with the numerical so- 
lutions for all reasonable parameter values. This expression highlights the 
connection between /jo„ and the thermal evolution of the Stromgren bub- 
bles: Essentially, the higher the quasi-static temperature, the lower the 
recombination rate A^ec, and the higher fion becomes. 

The value of fion relevant to computing the ionization redshift is obvi- 
ously that where z — ■Zjon- we have seen, Zjon typically lies between Izyif 
and Sz^jr- Substituting T 2, 500K into the expression for A^ec and taking 
Oo ~ 1 and z = Zion ~ "^-^Zmr: the above reduces to 



■""^ 1+0.8/1^6(1 + ^^^^)3/2' 

SO we see that fion will be of order unity unless Zyir » 15 or hQi, » 0.03. 
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discussions on the subject of the paper and W. Hu, A. Reisenegger, D. 
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Figure 1: Virialization redshifts for objects of varioTis masses. 
The virialization redshift, the redsliift at which the bulk of the objects of 
mass Mc form, is plotted for a number of cosmological models. In all cases 
shown, rio + Ao = 1 and Sc = 1.69. 
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Figure 2: Volume fraction ionized for various scenarios. 
The volume fraction of the universe that is in ionized Stromgren bubbles 
is plotted as a function of redshift for various parameter choices, corre- 
sponding to n = 1 CDM (optimistic, less optimistic, middle-of-the-road and 
pessimistic cases) and the tilted power spectrum (n=0.7) variant of CDM. 
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Figure 3: Ionization redshift for various scenarios. 
The redshift at which x = 0.5 plotted as a function of the net efficiency. 
The four curves correspond to four of the choices of Zyir in Table |l|: 41.4, 
17.2, 8.4 and 4.8 from top to bottom. 
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Figure 4: Opacity for completely ionized IGM. 
The Thomson opacity Ps{z). the probability that a CBR photon has been 
scattered at least once after the redshift z, is plotted for four different choices 
of /tOb for the case where the IGM is completely ionized at all times. For 
more realistic scenarios where ionization occurs around some redshift -Zjoru 
the opacity curves simply level out and stay constant for z 3> Zjon- 
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Figure 5: Last-scattering surface for completely ionized IGM. 
The probability distribution for the redshift at which a CBR photon was 
last scattered, the so called visibility function, is plotted for four different 
choices of hQh for the case where the IGM is completely ionized at all times. 
For more realistic scenarios where ionization occurs around some redshift 
Zion, the curves are unaffected for z <C Zion, vanish for Zion <^ z <^ 10^ and 
have a second bump around z ^ 10^. 
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Figure 6: Opacity for completely ionized IGM as function of angle. 
The total Thomson opacity Pg, the probability that a CBR photon has been 
scattered at least once since the recombination epoch, is plotted as a function 
of the angle in the sky that the horizon subtended at the reionization epoch. 
This is the largest angular scale on which fluctuations can be suppressed. 
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Figure 7: Last-scattering surface for completely ionized IGM as function of 
angle. 

The probability distribution for the angle subtended by the horizon when a 
CBR photon was last scattered, the so angular visibility function, is plotted 
for four different choices of hQ^ for the case where the IGM is completely 
ionized at all times. For more realistic scenarios where ionization occurs 
around some redshift Zion, corresponding to an angle Oiom the curves are 
unaffected for » 9ion, vanish for 2° <^ 9 <^ 9 ion and have a second bump 
around ~ 2°, the horizon angle at recombination. 
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Figure 8: Total opacity for various models. 
The total opacity, the probability that a CBR photon has been scattered at 
least once since the recombination epoch, is plotted for a variety of models. 
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Figure 9: Temperature evolution in intergalactic Stromgren bubbles. 
The temperature evolution is plotted for IGM exposed to a UV flux strong 
enough to keep it completely photoionized. In this example, h = 0.5, 
Qfe = 0.06, and T* = 36,900. The upper dashed line is T* , the temper- 
ature corresponding to the average energy of the released photoelectrons, 
towards which the plasma is driven by recombinations followed by new pho- 
toionizations. The lower dashed line is the temperature of the CBR photons, 
towards which the plasma is driven Compton cooling. The three solid curves 
from left to right correspond to three different redshifts for becoming part 
of a Stromgren bubble. The first time the hydrogen becomes ionized, its 
temperature rises impulsively to T*/2. After this, Compton cooling rapidly 
pushes the temperature down to a quasi-cquilibrium level, where the Comp- 
ton cooling rate equals the recombination heating rate. 
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Figure 10: Ionization efficiencies for various scenarios. 
The ionization efficiency, the fraction of the UV photons that produce a 
net ionization, is plotted for three different parameter combinations. In 
all cases, T* = 36, 900K, the value appropriate for the radiation from the 
population 3 star in Table |^. The solid lines are the exact results from 
numerical integration of equation (^). The dotted lines are the analytic 
fits, which are seen to agree well in the redshift range of interest, which is 
typically z twice or three times z^ir- 
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